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Motivated by results of recent analytic studies, we present 
a numerical investigation of the late-time dynamics of scalar 
test fields on Kerr backgrounds. We pay particular attention 
to the issue of mixing of different multipoles and their fall-off 
behavior at late times. Confining ourselves to the special case 
of axisymmetric modes with equatorial symmetry, we show 
that, in agreement with the results of previous work, the late- 
time behavior is dominated by the lowest allowed Z-multipolc. 
However the numerical results imply that, in general, the late- 
time fall-off of the dominating multipole is different from that 
in the Schwarzschild case, and seems to be incompatible with 
a result of a recently published analytic study. 



The purpose of the studies presented in this paper is to 
point out the need for further work towards a definitive 
answer to the following question: How do perturbations 
of a rotating black hole get radiated away? 

Using the Teukolsky formalism (lj, we can search for 
an answer to that question by studying the evolution of 
a complex spin-weighted wavefunction that is defined in 
terms of curvature perturbations of a fixed black hole 
background geometry. The late-time dynamics of gravi- 
tational perturbations is governed by features very sim- 
ilar to those occurring for scalar test-fields (^||. In the 
special case of vanishing spin- weight, the Teukolsky equa- 
tion reduces to the linear wave equation for a scalar test 
field $ on a Kerr background, which has a somewhat 
simpler structure than the Teukolsky equation. For the 
sake of simplicity, the present discussion is confined to 
the evolution of scalar test fields. 

For a spherically symmetric, asymptotically flat back- 
ground geometry, the angular variables can be separated 
off, and with an appropriate radial coordinate x the wave 
equation for a scalar test field <!> reduces to a (1 + 1) 
dimensional wave equation with an effective potential 
V(x), 



[d tt -d xx + V(x)]*(t,x)=0 



(1) 



When observed from a fixed location x = R, the dynam- 
ics of the wave consists of three stages (cf. Fig. |l|): During 
the initial burst phase (i), the waveform is dominated by 
the initial data, in our case given by an outgoing pulse 
with compact support. The initial burst is followed by 
the exponentially damped quasinormal ringing (ii) and 
the so-called tail phase (Hi), in which the amplitude of 
the field $ falls off according to a power-law, |$| cx 





(b) 



FIG. 1. Wave propagation on a curved background: (a) 
Three stages of the dynamics: initial burst phase (i), quasi- 
normal ringing (ii) and tail phase (Hi), (b) The quasinormal 
ringing and the tail phase are understood in terms of scatter- 
ing and interference - indicated by the dotted lines - of the 
wave off an effective potential V(x), depicted by the dashed 



The discovery that waves propagating on a curved 
background die off in tails goes back to the original work 
by R. Price in the early 1970s M. In the special case 
of a background geometry given by a nonrotating black 
hole, and when the field is not initially static, the late 
time decay at a fixed location is given by |$| cx £~( 2 '+ 3 ) ; 
where I is an angular separation parameter. In the last 
three decades, power-law tails on spherically symmetric 
backgrounds have been studied extensively by both an- 
alytic and numerical methods not only at a fixed 
location (i.e. at timelike infinity), but also at the horizon 
of the black hole and at future null infinity. 

This situation is in dramatic contrast to the case of 
a background given by a rotating Kerr black hole, in 
which until very recently no analytic predictions for the 
tail behavior were available. In fact it had been quite 
unclear whether tails are present at all, and, if they are, 
what the value of the power-law coefficient is. Numer- 
ical work, however, was done in an attempt to general- 
ize the above results to the Kerr case ||. Due to the 
breakdown of spherical symmetry of the background ge- 
ometry, a multipole decomposition in terms of spherical 
harmonics loses its special meaning and it is no longer 
possible to expand <I> in spherical harmonics in the usual 
way. Although the functions Y™ still form a complete 
set of functions for the description of a generic angular 
dependence, one now has to take into account the inter- 
action between different Z-multipoles in the course of the 
time evolution. Modes with different values of to, how- 
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ever, are not coupled because the background geometry 
is axisymmetric. Thus in order to take into account the 
contributions of all the different multipoles occurring in 
the numerical time evolution, only the azimuthal angular 
coordinate is separated off, and one has to solve a (2 + 1) 
dimensional evolution problem for the field 'J, defined by 

$ = *(t,r*,0)e im ^ , (2) 

where r* is the Kerr tortoise coordinate and is the 
azimuthal Kerr coordinate, that was chosen for technical 
reasons described in |§) . 

The main conclusion from the examples described in 
H can be stated in the following way: For a fixed value 
of the azimuthal separation parameter m, the late time 
dynamics is dominated by the lowest allowed ^-multipole 
that is compatible with the choice of m (I has to satisfy 
the condition I > m) and with the equatorial symmetry 
properties of the initial data. Note that symmetric and 
antisymmetric modes do not get mixed during the evolu- 
tion. The fall-off behavior of the dominating multipole is 
given by |$| oc £~( 2 '+ 3 ). This statement is illustrated in 
Fig. ^| by a specific example in which the initial pulse is 
given by I — 2, m = 0, which is not the lowest axisymmet- 
ric mode that is symmetric with respect to the equatorial 
plane: During the evolution there is a transition to the 
I — m — mode. The numerically determined power-law 
coefficient governing the decay of the I = m = mode is 
fi = 3.15 ||. What, in this logarithmic diagram, appears 
to be a notch for angles smaller than 9 ps 0.9, in fact rep- 
resents the change of sign that the field $ experiences for 
those angles during the transition from the initial I = 2 
distribution to the / = mode. 




FIG. 2. Plots of |$| for m = 0, a = 0.9999 for different 
angular directions < 6 < tt/2. The initial pulse was given 
by I — 2. The time t and |$| are shown on a logarithmic 
scale. For late times the numerically determined power-law 
coefficient is fj, = 3.15. 

The example shown in Fig. ^| illustrates the late-time 
behavior for one special choice of parameters. However 
there are several questions that remain open: What hap- 
pens, for example, if we take initial data given by an 



I = 4 multipole? In analogy to the situation depicted in 
Fig. we expect a transition to the I = mode that will 
decay as t~ 3 . More generally, is seems obvious that any 
axisymmetric initial pulse with equatorial symmetry will 
descend to I — and that the late-time power-law decay 
of the field <£> will be given by |$| oc t~ 3 . 

Recently, these and other issues have been addressed 
in two different analytic studies of late-time phenomena 
occurring in the course of the evolution of a massless 
scalar field on rotating black hole backgrounds. 

L. Barack and A. Ori [jl0| have presented an analysis 
that is done in the time domain and yields expressions 
for the asymptotic behavior at future null infinity, at the 
black hole horizon, and at timelike infinity. After decom- 
posing the field <I> into spherical harmonics, 

$ = (r 2 + a 2 )- 1/2 E *'"*(*> r ) y "^' ' 

I'm 

they consider an expansion of W m (t, r) in inverse powers 
of advanced time (late-time expansion). According to 
their prediction, to leading order in t, the late-time decay 
of a mode specified by I and m at timelike infinity is given 
by 

$ K t-i-M-3-9 , (3) 

where q = if I + m is even, and 9 = 1 otherwise. Note 
that this prediction by Barack and Ori applies only in 
the case in which the I = mode is present in the initial 
data. 

Independently, S. Hod has carried out an analysis that 
is based on a frequency decomposition of the scalar field 
$ pT[ . His analysis leads to the following conclusion 
for the late-time behavior at timelike infinity: The late- 
time behavior at timelike infinity is dominated by the 
lowest allowed multipole, I — \m\, if I* — \m\ is even, 
and / = \m\ + 1 otherwise, where I* denotes the initial 
multipole. The late-time fall-off is given by 

$oct _, *- |m| - p - 1 ifT>H + 2, (4a) 
$oci -2r-3 if I* = jm|, |m| + 1 , (4b) 

where p = if I* — \m\ is even, and p = 1 otherwise. 

For certain special choices of parameters, the expres- 
sions ([|) and (||) are in agreement. Consider for example 
the case of an initial pulse given by I — m = 0. Barack 
and Ori, as well as Hod, predict a late-time fall-off given 
by |$| oc t~ 3 . In this case both analytic predictions also 
agree with the numerical results, shown in [||. 

In general, however, a comparison of the expressions 
(||) and ((J) is made difficult by the fact that, as pointed 
out above, Eq. (||) does not allow a general statement 
about the dominating late-time behavior. In contrast to 
(^), Hod's formula ([|) makes a definite statement about 
the dominating late time decay, depending on the initial 
multipole. If one, for example, sets m = and chooses an 
initial pulse with /* = 4, Hod's analysis implies |<f>| oc t~ 5 . 
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According to Hod's analysis, the I — multipole "remem- 
bers" the initial configuration I* in the sense that its de- 
cay rate does not depend only on the choice of m and 
the equatorial symmetry properties of the initial data, 
but also on I*. 

The presence of two different analytic studies, one of 
which contains a rather surprising prediction, has moti- 
vated us to carry out a numerical investigation of the late 
time fall-off behavior. 

For the axisymmetric evolution studies presented in 
this paper, we use the numerical (2 + 1) evolution code 
from H with the black hole mass set to unity. We made 
several modifications to the previously used code. In- 
stead of using the Boyer-Lindquist 9 coordinate, we are 
using y = — cos 9, because this coordinate choice allows 
us to keep track of the behavior of the various angu- 
lar modes more accurately. The y derivatives of \& are 
approximated using fourth order finite differences. We 
have used quadruple precision floating point arithmetic 
for the numerical computations, so that we can follow the 
dynamics up to late evolution times. This issue is not of 
importance for low initial multipoles, but for I = 4 and 
higher it becomes relevant because of the steep fall-off 
at the beginning of the tail regime. The accuracy of the 
code was tested using standard convergence tests, based 
on numerical results obtained for short evolution times. 

For m — 0, in the t,r*,y coordinates, the wave equa- 
tion for \I>(t, r*,y) reads 



2rA 



(r 2 + a 2 f 



S 2 



A(l-y 2 ) 

£ 2 



(5) 



where we are using the standard abbreviations A = r 2 — 
2Mr + a 2 , and £ 2 = (r 2 + a 2 ) 2 - a 2 A(l -y 2 ) . The initial 
data are given by the Gaussian-like, initially outgoing 
pulse with compact support given by 50M < r* < 150M 
used in g. 

The following figures represent an analysis of the data 
obtained from numerical evolutions. For a range of time 
parameter values t, the signal VP at fixed radius is decom- 
posed into angular modes. We start the evolutions with 
modes that are symmetric with respect to the equato- 
rial plane (7 = or 4). Since symmetric modes do not 
get mixed with antisymmetric modes, we can discard the 
latter and just decompose the signal into even modes. 
Since we expect the signal at late times to be dominated 
by the lowest multipoles, we only take into account the 
three lowest even modes, 



ty(t,r* = const., y) = 
a(t)P (y) + p(t)P 2 (y) 



(6) 



where the i"b,2,4(j/) are the standard Legendre polynomi- 
als. 

In Fig. H we depict the fall-off behavior of the absolute 
values of the expansion coefficients \a\, \/3\, and I7I, ex- 
tracted at r* = 10, as functions of the Boyer-Lindquist 



time t for the test case a = 0. The grid spacings are given 
by dr* = 0.04 and dy = 0.01. The initial pulse is given 
by the I = 4 multipole, and the computation is done with 
the same (2 + 1) evolution code that is used in the a^0 
case. Since for a = the background is spherically sym- 
metric, the I = 4 multipole has to be preserved during 
the evolution. Hence we know that all contributions from 
modes other than I = 4 are numerical artifacts. A least 
squares fit of a function f(t) = ci _AI to j(t) yields the 



power-law j(t) = -3.33xl0 18 T 



The power-law co- 



efficient [i is within less than 0.7% in agreement with the 
theoretical value of fi = 11. The expansion coefficients of 
the numerical contamination modes I = 0, 2 are orders of 
magnitude smaller than the I = 4 mode, hence establish- 
ing the accuracy of the numerical evolution method. 
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FIG. 3. Log-log plots of the coefficients a, (i, and 7 for 
m = 0, a = 0, evaluated at r* = 10. The initial data are 
given by an I = 4 pulse. For late times, the decay of the I = 4 
mode is governed by the power-law j(t) cx t~ 10 93 . 

What can we expect for a ^ 0? Intuitively we might 
anticipate a transition to the I = mode that will decay 
as t~ 3 , similar to the situation illustrated in Fig. ||. How- 
ever we know that Hod's analytic prediction, which we 
briefly discussed above, yields | \l/ 1 oc i~ 5 for the dominat- 
ing I = mode at late times, contrary to a naive guess. 
For grid spacings given by dr* = 0.04 and dy = 0.01, a 
fit to the numerically computed a(t) yields the power-law 
a(t) oc i~ 6 - 04 for the decay of the I = mode which is 
dominating at late times. This result is deeply surprising. 
The transition to the I = mode is quite intuitive, but we 
might have expected a decay rate of i -3 , independently 
of the particular choice of the initial value of I, as long as 
the initial pulse is symmetric with respect to the equa- 
tor. The numerical result is not only counter-intuitive, 
it also seems to be incompatible with the analytic pre- 
diction given in pj] ]. This suggests that the computed 
a(t) oc t~ e ' 04 behavior might represent a numerical arti- 
fact rather than a physically meaningful result. 
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FIG. 4. Plots of fi a for m = 0, a = 0.9999, observed at 
r* = 10 (boxes) and r* = 25 (stars). The initial pulse was 
given by I = 4. 

This suspicion is confirmed by Fig. |I| that shows the 
grid size dependence of the numerically extracted power- 
law coefficient /i Q , that governs the decay of a at r* = 10 
(boxes) and at r* = 25 (stars) The results were ob- 
tained from four runs with dr* = 0.16,0.1,0.08,0.064 
and dy = 0.01,0.00625,0.005,0.004, i.e. the discretiza- 
tion in the y direction was chosen to be finer than in the 
cases above. To limit the computational expenses, the 
run times were shorter than in the previously considered 
cases, and fi a was extracted for 600 < t < 800. As illus- 
trated in Fig. ^, a fit of a function fi a = c\ + ci dy C3 to 
the numerical data yields a value of C3 that is close to 4, 
indicating that the discretization error is dominated by 
the finite difference approximation of the y derivatives. 
The continuum limit of fi a , however, seems to be neither 
5 nor 6 and shows a rather weak dependence on the loca- 
tion of the observer, possibly a numerical artifact. Our 
numerical work strongly suggests that an extension of the 
existing analytic and numerical studies is necessary. 

The studies presented in this paper show that, for m 
set to zero and initial data with equatorial symmetry, 
the late-time behavior of scalar test fields on Kerr back- 
grounds is dominated by the lowest allowed Z-multipole, 
namely 1 = 0. The numerical results imply that, in gen- 
eral, the late-time fall-off of that dominating multipole 
is different from that in the Schwarzschild case, and ap- 
pears not to agree with a central result of Hod's recent 
analytic study [fj"lf . Compatibility of our numerical re- 
sults with the analysis of Ori and Barack |Ic| |, requires 
an extension of their studies to initial data that do not 
contain the I = multipole. The characterization that 
the initial pulse is "pure I — 4" is ineluctably related to 
the use of the Boyer-Lindquist 9 coordinate. Without 
explicit reference to the Boyer-Lindquist coordinate sys- 
tem, the numerical result can be stated as: There is an 
axisymmetric initial pulse with equatorial symmetry that 
exhibits a transition to I — and the late-time power-law 
decay of the field $ is given by |$| cx t~^, where the nu- 



merically determined power-law coefficient fi is different 
from the Schwarzschild value /1 = 3. Additionally, our re- 
cent studies suggest that conclusions about the late-time 
power-law decay, made in [pi on the basis of particular 
examples, are not valid for generic situations. 

The numerical results suggest that further analytic 
work is needed. It will also be useful to perform numer- 
ical computations up to later evolution times to under- 
score the validity of the statements made in this paper. 
Future numerical work should include studies of higher, 
non-axisymmctric multipoles, and also the physically rel- 
evant case of gravitational perturbations. 
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